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Abstract. Excitable cellular automata with dynamical excitation interval exhibit a 
wide range of space-time dynamics based on an interplay between propagating excita- 
tion patterns which modify excitability of the automaton cells. Such interactions leads 
to formation of standing domains of excitation, stationary waves and localised excita- 

^N) tions. We analysed morphological and generative diversities of the functions studied 

f-^ and characterised the functions with highest values of the diversities. Amongst other 

^ intriguing discoveries we found that upper boundary of excitation interval more signifi- 

cantly affects morphological diversity of configurations generated than lower boundary 
of the interval does and there is no match between functions which produce configura- 

^ ^ tions of excitation with highest morphological diversity and configurations of interval 

boundaries with highest morphological diversity. Potential directions of future studies 

I— I of excitable media with dynamically changing excitability may focus on relations of 

the automaton model with living excitable media, e.g. neural tissue and muscles, novel 

r \ materials with memristive properties, and networks of conductive polymers. 

Keywords: excitation, automata, diversity, localisations, patter formation 

' 1. Introduction 

Since their popularisation in [5j, excitable cellular automata became a convenient tool 
QQ for studying complex phenomena of excitation dynamics and chemical reaction-diffusion 

\^ activities in physical, chemical and biological systems [7, 3J. The cellular automata offers 

quick 'prototyping' of complex spatially extended non-linear media. The examples of 
'best practice' include models of Belousov-Zhabotinsky reactions and other excitable 
^ systems [H [H], chemical systems exhibiting Turing patterns [121 El HD], precipitating 

systems |2], calcium wave dynamics [11], and chemical turbulence [6]. 
I In a classical Greenberg-Hasting [5| automaton model of excitation a cell takes three 

states — reseting, excited and refractory. A resting cell becomes excited if number 
of excited neighbours exceeds a certain threshold, an excited cell becomes refractory, 
rN and a refractory cell returns to its original resting state. In 1998 p!], we introduced an 

excitable cellular automaton, where a resting cell is excited if a number of its excited 
neighbours belongs to some fixed interval [^1,^2]- The interval [^1,^2] was called an 
excitation interval. For a two-dimensional cellular automaton with eight-cell neighbour- 
hood boundaries of the excitation interval satisfy the condition: 1 < ^1 < ^2 ^ 8. We 
found that by tuning 9i and 62 we can force the automaton to imitate almost all kinds 
of excitation dynamics, from classical target and spiral waves observed in physical and 
chemical excitable media to wave-fragments inhabiting sub-excitable media. 
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How does excitation dynamics change if we allow boundaries of the excitation interval 
to change during the automaton development? We partially answer the question in 
present paper by making the interval [^i(x), 02(^)] of every cell x to be dynamically 
updatable at every step t depending on state of the cell x and numbers of excited and 
refractory neighbours in the cell x's neighbourhood. 

The excitable automata with dynamical excitation intervals are defined in Sect. [2} 
Morphological diversity of cellular automata (measured using Shannon entropy and 
Simpson index) with different functions of interval updates is characterised in Sect, [sj 
Section |4] characterises generative diversity (measured in terms of different configurations 
generated during space-time development of automaton starting with a single non-resting 
cell) of the local transitions. Some afterthoughts are offered in Sect, p] 



2. Dynamical excitation intervals 

Let X* and x''~^^ be states of a cell x at time steps t and t + 1, and cr?j.(x) be a sum of 
excited neighbours in cell x's neighbourhood u{x) = {y ■.\x — yji^ = 1}. Cell x updates 
its state by the following rule: 

'+, ifx* 



• and a\{x)+ S [e\{x) , ei{x)\ 



if X* 



+ 



otherwise 



A resting cell is excited if number of its neighbours belongs to excitation interval [^i(x), 0\{x)\, 
where 1 < 6\{x)^9\{x) < 8. The boundaries 9\{x) and 9\{x) are dynamically updated 
depending on cell x's state and numbers of x's excited o"+(x) and refractory o"L(x) neigh- 
bours. A natural way to update boundaries is by increasing or decreasing their values 
as follows: 

e\+\x) = i{e\{x) + ^, 
= e(^*(x) + A2^ 



(^))) 



where 



Ai 



Ti , if X 
Ta, ifx 
0, if X = 



+ 
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T2, 



if X 
if X 



+ 



0, if X = 



and (j){a — 6) = 1 if a > 6, if a = 6 and -1 if a < 6, and ^(a) = 1 if a < 1 and 8 if a > 8. 
Boundaries of excitation interval [6\{x)^ ^K^)] ^^^^ updated independently of each other. 
Local excitation rules are determined by values of Ti,- ■ - Ta. We therefore address the 
functions as tuples E{TiT2T^T4) which range from £'(—1 — 1 — 1 — 1) to £'(1111). 

Functions £(0060), a G { — 1,0,1} represent rules with fixed upper boundary 62 of 
excitation and dynamically updated lower boundary ^1. These are equivalent to dy- 
namically updated thresholds of excitation. Functions £(OaOa), a G { — 1,0, 1} represent 
rules with fixed lower boundary and dynamical upper boundary of excitation interval. 

The experiments are conducted on a cellular array of n x n cells with periodic boundary 
conditions. In a typical experiment we perturb resting cellular array with a localised 
domain of excitation, wait till transient period is over (1000 iterations is enough) and 
most excitation patterns collide and disappear and persist indefinitely, and then analyse 
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three configurations: configuration of excitation represented by an array of cells states 
X*, and configurations of lower 9i{x) and upper 02ix) boundaries of excitation intervals. 

Initially ^2(^) = ^ for any x. In experiments we considered initial conditions 9i{x) = 1 
and 0^{x) = 2. The following scenaria of initial excitation were played: 

• (++)-start, 9i = 2: all cells are resting but two neighbouring cells are assigned 
excited state, 

• Rl-start: let D be a disc radius n/4 centred in the array L oinxn cells, all cells 
are resting but cells lying in D are assigned excited states with probability 0.2 
and 9\{x) = 1 for any x, 

• R2-start: all cells are resting but cells lying in D get excited states with proba- 
bility 0.2 and 9\{x) = 2 for any x, 

• Dl-start: all cells are resting but cells lying in D get excited states with proba- 
bility 0.1 or refractory states with probability 0.1 and 9\{x) = 1 for any x, 

• D2-start: all cells are resting but cells lying in D get excited states with proba- 
bility 0.1 or refractory states with probability 0.1 and 9\{x) = 2 for any x, 

• (— +)-start: all cells are resting but one cell is excited and its western neighbour 
is refractory, 

• ( — h +)-start all cells are resting but one cell is excited, its first order western 
neighbour is excited and its second order western neighbour is refractory. 

Cell states were represented by colours and grey levels as follows: excited state + is 
red (c. 76 grey), resting state is white and refractory state — is blue (c. 28 grey). Colour 
values of excitation interval boundaries 9i and 92 are following: 1 is white, 2 is green or 
150 grey, 3 is yellow or 226 grey, 4 is blue or 28 grey, 5 is magenta or 104 grey, 6 is cyan 
or 178 grey, 7 is red or 76 grey, and 8 is black. 

3. Morphological diversity 

We evaluated morphological diversity of configurations of excitation and using Shan- 
non entropy and Simpson's index. Let W = {o, +, — } be a set of all possible con- 
figurations of a 9-cell neighbourhood w{x) = u{x) U x, x G L. Let c be a con- 
figuration of automaton, we calculate number of non-resting neighbourhood configu- 
rations as r/ = X^2:gL^(^)' ■^here e(x) = if for every resting x all its neighbours 
are resting, and e(x) = 1 otherwise. The Shannon entropy is calculated as f = 
~ S«)eH/(^(^)/^ • ln{u{w) /rf)), where vivo) is a number of times the neighbourhood 
configuration w is found in automaton configuration c. Simpson's index is calculated as 
5 = 1 — X]«;eVK(^(^)/^)^- '^^^ measures £ and S were calculated on configurations of 
cell-states and interval boundaries after long transient period, sufficient enough for any 
perturbation to settle down. 

The diversity of excitation patters is evaluated using S-£ plots. See examples of S-£ 
plots for Dl-start, 6\{x) = 1, in Fig.[l]and (++)-start, 6\{x) = 2, in Fig.[2| Distributions 
of functions by their values of Shannon entropies for 9i and 92 are illustrated in Figs. [3] 
andH 
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Figure 1. Morphological diversity of functions for Dl-start, 6i{x) = 1 
for all x: Simpson's index S (horizontal axis) vs Shannon entropy £ (ver- 
tical axis) for configuration of excitable array of 200 x 200 cells with peri- 
odic boundary condition, recorded at t = 1000. Encircled data points are 
seven functions with highest morphological diversity specified in column 
Dl-start, S-£ in Tab. [3j 
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Figure 2. Morphological diversity of functions for (++)-start, 6i{x) = 2 
for all x: Simpson's index S (horizontal axis) vs Shannon entropy £ (ver- 
tical axis) for configuration of excitable array of 200 x 200 cells with peri- 
odic boundary condition, recorded at t = 1000. Encircled data points are 
seven functions with highest morphological diversity specified in column 
(+-|-)-start, S-£ in Tab.jsj 
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Figure 3. Morphological diversity of functions for Dl-start, 6i{x) = 1 
for all x: Shannon entropy £1 for configuration of 61 (horizontal axis) 
vs Shannon entropy £2 for configuration of 82 (vertical axis), recorded 
at t = 1000. Encircled data points are seven functions with highest 
morphological diversity specified in column Dl-start, £1-82 in Tab. [3) 
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Shannon entropy of 9i vs Shannon entropy of 02 
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Figure 4. Morphological diversity of functions for (++)-start, 6i{x) = 2 
for all x: Shannon entropy £i for configuration of 6i (horizontal axis) 
vs Shannon entropy £2 for configuration of O2 (vertical axis), recorded 
at t = 1000. Encircled data points are seven functions with highest 
morphological diversity specified in column (++)-start, £1-82 in Tab. [sj 
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Top seven functions showing highest values of diversity indices, e.g. those encircled 
in examples Fig. [l]-[4j are grouped in Tab. |3]for various scenarios of initial start. 

Finding 1. There is no match between functions which produce configurations of excita- 
tion with highest morphological diversity and configurations of interval boundaries with 
highest morphological diversity.. 

Amongst functions listed in Tab. [s] only function E{\ — 101) gets into top seven 
functions with highest diversity of both excitation and interval boundaries for scenario 
( — h +)-start. Exemplar configurations of excitation and interval boundaries generated 
by automaton governed by E{\ — 101) are shown in Fig. [To| 3. Function E[l — 101) is 
also amongst functions with highest diversity for D2- and ( — |-)-start. The function 
governs the following update of the excitation interval boundaries. Low boundary Oi{x) 
of excitation of cell x is updated only cell x is excited. The boundary 9i{x) increases if 
cell X has more excited neighbours than refractory neighbours, 6i{x) decreases if number 
of refractory neighbours of x exceeds number of excited neighbours. The boundary Oi{x) 
does not change if cell x has the same number of excited neighbours as refractory neigh- 
bours. Upper boundary 02{x) increases if cell x's dissents with excitation-refractoriness 
ratio in its neighbourhood: x* = + and a'\^{x) < a^_{x) or = — and (7^^{x) > (t'l(x). 
The boundary 92{x) decreases if cell x conforms to excitation-refractoriness ratio in its 
neighbourhood: x* = + and cr^{x) > a^_{x) or x* = — and cr^{x) < (t^_{x). Increase 
of 9i and decrease of 62 lead to decrease cell's excitability. Thus we can characterise 
function E{\ — 101) as follows: excitability of a cell decreases if the cell dissents with its 
neighbourhood and increases otherwise. 

Finding 2. Functions E{—\ — l{)l) and E{ — 1{}1\) generate configurations of excitation 
interval boundaries with highest morphological diversity for three and two types of initial 
stimulation, respectively, and functions E (000 — 1) and E (0100) generate configurations 
of excitation with highest morphological diversity for two types of initial stimulation. 

The function £'(—1 — 101) generates configurations with highest morphological diver- 
sity of ^1 and 02 configurations for R2-, Dl- and (— +)-starts, see examples in Fig. [Tja and 
|9j3, and £^(—1011) generates highest morphological diversity configurations of 61 and 62 
for R2- and D2-starts (Figs.[6]o and[8)D). The functions £(000-1) and £(0100) produce 

highest morphological diversity configurations of excitation for (++)- and ( |-)-starts 

(£(000 - 1)) and R2- and D2-starts (£(0100)). See examples for £(000 - 1) in Figs.jsji 
and[T0}i and £(0100) in Figs. [6ji and [sji. 



Finding 3. Function £(0100) generates most morphologically diverse excitation patterns 
in larger, comparing to other functions, number of initial conditions. 

Function £(0100) is amongst top seven functions with highest morphological diversity 
of excitation in (++)-, R2-, D1-, D2- and ( — h +)-starts (Tab. |3]). The function gen- 
erates most morphologically diverse excitations in D2-start. Examples of configurations 
generated by £(0100) are shown in Figs. |6^ and [8^. The function £(0100) shows how 
dynamics of excitation can be tuned by changing only upper boundary of the excitation 
interval with lower boundary fixed. Value 9i{x) is not updated. Value 92{x) increases 
if cell X is excited and it has more excited than refractory neighbours, the value 92{x) 
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(a) ++-start, £;(000 - 1) 




(b) ++-start, E{-11 -1-1) 




(c) Rl-start, E{1 - 111) 



Figure 5. Examples of most morphologically diverse configurations gen- 
erated in (++)-start (ab) and Rl-start (c). (a) Configurations with high- 
est morphological diversity of excitation generated by function £'(000—1). 
(b) Configurations with highest morphological diversity of interval bound- 
aries 9i and 62 generated by function £(—11 — 1 — 1). (c) Configurations 
with highest morphological diversity of interval boundaries 61 and 62 
generated by function E{1 — 111). Configurations of 61 (left), excitation 
(middle) and 62 (right) are taken in 200 x 200 cell array, at t = 1000. 
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(b) R2-start, _B(-1011) 



Figure 6 . Examples of most morphologically diverse configurations gen- 
erated in R2-start. (a) Configurations with highest morphological diver- 
sity of excitation generated by function E{OWO). (b) Configurations with 
highest morphological diversity of interval boundaries 9i and 02 generated 
by function £^(—1011). Configurations of Oi (left), excitation (middle) 
and 02 (right) are taken in 200 x 200 cell array, at t = 1000. 

decreases if cell x is excited and has more refractory than excited neighbours. Excitabil- 
ity of a cell decreases if the cell dissents with excitation ratio in its neighbourhood, and 
increases otherwise. 

Finding 4. Function E{—10 — 1 — 1) generates most morphologically diverse patterns of 
interval boundaries in larger, comparing to other functions, number of initial conditions. 

Function E{—10 — 1 — 1) gets in top seven functions with highest morphological di- 
versity of 6*1 and 62 patterns in (++)-, R2-, D1-, D2-, and ( — |-)-starts (Tab. [s]). In 
automata, governed by this function, 62{x) is not updated if cell x is excited. Otherwise, 
^1(2;) and 62{x) increase if refractory neighbours outnumber in the cell x's neighbourhood 
and decrease if excited neighbours outnumber refractory neighbours. 

Finding 5. Upper boundary of excitation interval more significantly affects morpholog- 
ical diversity of configurations generated than lower boundary of the interval does. 
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(a) Dl-start, £(0101) 




(b) Dl-start, E{-1 - 101) 




(c) Dl-start, £(11 -1-1) 



Figure 7. Examples of most morphologically diverse configurations gen- 
erated in Dl-start. (a) Configurations with highest morphological diver- 
sity of excitation generated by function S(OIOI). (b) Configurations with 
highest morphological diversity of interval boundaries 9i and 62 gener- 
ated by function £^(—1 — 101). (c) Configurations with highest morpho- 
logical diversity of interval boundaries 9i and 62 generated by function 
£^(11 — 1 — 1). Configurations of 9i (left), excitation (middle) and 62 
(right) are taken in 200 x 200 ceh array, at t = 1000. 




(b) D2-start, £(-1011) 



Figure 8. Examples of most morphologically diverse configurations gen- 
erated in D2-start. (a) Configurations with highest morphological diver- 
sity of excitation generated by function E{OWO). (b) Configurations with 
highest morphological diversity of interval boundaries 9i and 02 generated 
by function £^(—1011). Configurations of Oi (left), excitation (middle) 
and 02 (right) are taken in 200 x 200 cell array, at t = 1000. 



There are only two functions, i?(1010) and i?(1000), where only lower boundary 0i is 
updated in Tab.[3j Function £'(1010) represents a situation when 0i{x) is independently 
of a state of cell x: 0i{x) increases if number of excited neighbours exceeds number of 
refractory neighbours, and 0i{x) decreases if refractory neighbours outnumber excited 
neighbours. In automata governed by function £^(1000) value of ^1(3;) is updated as 
above but only if cell x is excited. There are several functions with highest morphological 
diversity which represent fixed lower boundary and dynamical upper boundary, e.g. 
£;(000-l) (Figs. [§i and [Toji), £^(0100) (Figs. [6}i and [8}i), £;(0101) (Fig.[7|i), £(0-100), 
£;(0 - 101), and £(0001) (Tab.[3]). 



4. Generative diversity and localisations 



Generative diversity characterises how many different configurations are generated 
during space-time development of automaton starting with entirely resting configuration 
but single cell in a non-resting state. We consider two starting conditions: (++)-start 
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Figure 9. Examples of most morphologically diverse configurations gen- 
erated in (— +)-start. (a) Configurations with highest morphological di- 
versity of excitation generated by function £^(—1 — 100). (b) Configura- 
tions with highest morphological diversity of interval boundaries 6i and 
02 generated by function E{—1 — 101). (c) Configurations with highest 
morphological diversity of interval boundaries 9i and 82 generated by 
function £^(—1001). Configurations of 9i (left), excitation (middle) and 
62 (right) are taken in 200 x 200 cell array, at t = 1000. 
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(a) (- + +)-start, _E(000 - 1) 






(b) (- + +)-start, £(1 - 101) 





(c) (- + +)-start, i;(110 - 1) 



Figure 10. Examples of most morphologically diverse configurations 
generated in ( — |-+)-start. (a) Configurations with highest morphological 
diversity of excitation generated by function £'(000 — 1). (b) Configura- 
tions with highest morphological diversity of interval boundaries 6i and 
02 generated by function E(l — 101). (c) Configurations with highest 
morphological diversity of interval boundaries 9i and O2 generated by 
function £^(110 — 1). Configurations of Oi (left), excitation (middle) and 
02 (right) are taken in 200 x 200 cell array, at t = 1000. 
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Space filling vs Shannon entropy 
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Figure 11. Generative diversity of functions for (++)-start, O^^x) = 2 
for all x: Space filling (horizontal axis) vs Shannon entropy (vertical 
axis) for configuration of excitable array of 200 x 200 cells with peri- 
odic boundary condition, recorded at t = 1000. Encircled data points 
are seven functions with highest generative diversity specified in column 
(++)-start, in Tab. [2} 



and 0^{x) = 2 for any x; ( — |-)-start and O^^x) = 1; and, ( |-)-start and 9i{x) = 2. 

Generative diversity is evaluated using Shannon entropy and space-filling (a ratio of cells 
in a non-resting state). Functions generating configurations with maximum Shannon 
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(++)- and (- + +)-start, 6*? = 2 


(-+)-start, 6^? = 1 


-11-1-1 


-1 - 101 


01-1-1 


0-1-1-1 


100 - 1 


00-1-1 


1010 


01-1-1 


11 - 10 


01 - 11 


1110 


10-1-1 


1111 


11 - 11 



Table 2. Seven functions with highest generative diversity for (++)- 
and ( — h +)-start, Oi = 2 (first column) and ( — |-)-start, 9i = 1 (second 
column) . 



(++)-start ( — h)-start 



10- 1-1 

11- 1-1 

II - 10 
100-1 
101 - 1 
110-1 

III - 1 

1011 
1110 

Table 3. Functions supporting localizations in (++)- and ( — |-+)-start 
(left column) and ( — |-)-start (right column). 



-10 - 1 - 
-11 - 10 

00- 1-1 

01- 1-1 
01 - 10 

10- 1-1 

11- 1-1 
11 - 10 



entropy and minimum space-filling are assumed to have higher generative complexity, 



see example in Fig. 12 



Seven functions with highest generative diversity are listed in Tab. [2j Configurations 
of excitation and interval boundaries for £'(—11 — 1 — 1), ++-start, shown in Fig. [sja, 
E{—1 — 101), ( — |-)-start, in Fig.^ and configurations generated by functions -E(llll), 



E{ll - 10) and £(0 - 1 - 1 - 1) in Fig. 12 



Finding 6. Most localizations generated in (++)- and ( — h +)-starts are stationary. 

Around half of the functions generate configurations with localizations in case of R2- 
start (43), Dl-start (41), D2-start (43). We concentrate on functions which produce 
localizations in singleton starts. There nine functions for (++)- and eight function 
for ( — |-)-start , 6^ = 2 (Tab. [s]). Examples of configurations generated by functions 



£(00 - 1 - 1), £(101 - 1), £(1110) and £(00 - 1 - 1) are shown in Fig. 13 
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Figure 12. Examples of configurations generated by functions with 
highest generative diversity, (a) ^^(1111), (b) ^(11 - 10), (c) ^(0 - 
1 — 1 — 1). Automaton array has 200 x 200 cells, configurations of 6i 
(left), excitation (middle) and 02 (right) at t = 1000. 



Finding 7. Functions £^(11 — 10) and 100 — 1 are amongst top seven functions with high- 
est generative diversity supporting localised excitation dynamics in scenarios of (++)- 
start. Functions E{00 - 1 - 1), E{01 - 1 - 1), ^(10 - 1 - 1) and E{11 - 1 - 1) are 
amongst top seven functions with highest generative diversity supporting localised exci- 
tation dynamics in scenarios of (—+)-start 
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(d) (-+)-start, = 1, £(00 -1-1) 



Figure 13. Examples of configurations with localised excitations devel- 
oped in (++)-start (abc) and ( — |-)-start (d) scenarios. Size of cellular 
array is 200 x 200 cells, configurations of 6i (left), excitation (middle) 
and 02 (right) at t = 1000. 
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See Tabs. [2] and [3j Configurations generated by function E{11 — 10) are exemplified in 
Fig. [12)3 and function E{00 - 1 - 1) in Fig. [Tsji. Functions E{WO - 1) and ^(00 - 1 - 1) 
are the functions with minimal updates of excitation interval boundaries. In automata 
governed by function E{00 — 1 — 1) Oi{x) and 62{x) are updated only if the cell x 
is in refractory state: both boundaries decrease if excited neighbours of x outnumber 
refractory neighbours, and they increase refractory neighbours dominate. In automata 
governed by function £^(100 — 1) the boundary 9i{x) increases if the excited cell x has 
more excited neighbours than refractory ones, and boundary and decreases if the cell has 
more refractory neighbours. The boundary 92{x) is updated only if cell x is refractory: 
02{x) decreases if excited neighbours outnumber refractory neighbours, and it increases 
otherwise. 



5. Summary 

Excitable cellular automata with dynamical excitation interval exhibit a wide range 
of space-time dynamics based on an interplay between propagating excitation patterns 
which modify excitability of the automaton cells. Such interactions leads to formation of 
standing domains of excitation, stationary waves and localised excitations. We analysed 
morphological and generative diversities of the functions studied and characterised the 
functions with highest values of the diversities. Amongst other intriguing discoveries we 
found that upper boundary of excitation interval more significantly affects morpholog- 
ical diversity of configurations generated than lower boundary of the interval does and 
there is no match between functions which produce configurations of excitation with 
highest morphological diversity and configurations of interval boundaries with highest 
morphological diversity. Potential directions of futures studies of excitable media with 
dynamically changing excitability may focus on relations of the automaton model with 
living excitable media, e.g. neural tissue and muscles, and novel materials with memris- 
tive properties, and networks of conductive polymers. 
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